home *** CD-ROM | disk | FTP | other *** search
/ CU Amiga Super CD-ROM 19 / CU Amiga Magazine's Super CD-ROM 19 (1998)(EMAP Images)(GB)[!][issue 1998-02].iso / CUCD / Programming / LEDA / incl / LEDA.020+881 / Float.h < prev    next >
C/C++ Source or Header  |  1994-08-05  |  2KB  |  96 lines
  1. /*******************************************************************************
  2. +
  3. +  LEDA  3.1c
  4. +
  5. +
  6. +  Float.h
  7. +
  8. +
  9. +  Copyright (c) 1994  by  Max-Planck-Institut fuer Informatik
  10. +  Im Stadtwald, 6600 Saarbruecken, FRG     
  11. +  All rights reserved.
  12. + 
  13. *******************************************************************************/
  14.  
  15.  
  16. #include <LEDA/basic.h>
  17. #include <LEDA/Int.h>
  18. #include <math.h>
  19.  
  20. const double eps0 = ldexp(1,-53);  // 2^-53
  21.  
  22.  
  23. enum { ZERO = 0, NON_ZERO = 1, NO_IDEA = 2 };
  24.  
  25.  
  26. class Float {
  27.  
  28.   double num;
  29.   double mes;
  30.  
  31. public:
  32.  
  33. Float(Int i) 
  34. { if (i == 0) 
  35.      { num = 0; 
  36.        mes = 0; 
  37.       }
  38.   else 
  39.      { num = Itodouble(i);
  40.        int exp = 1 + lg(i); //Ilog(abs(i) - 1);
  41.        mes = ldexp(1, exp );
  42.       }
  43. }
  44.  
  45. Float(double d, double m) { num = d; mes = m; }
  46.  
  47. Float() { num = 0; mes = 0; }
  48.  
  49.  
  50. operator double() const { return num; }
  51.  
  52.  
  53. friend Float operator+(const Float& a, const Float& b)
  54. { return Float(a.num + b.num, 2*((a.mes > b.mes) ? a.mes : b.mes)); }
  55.  
  56. friend Float operator-(const Float& a, const Float& b)
  57. { return Float(a.num - b.num, 2*((a.mes > b.mes) ? a.mes : b.mes)); }
  58.  
  59. friend Float operator*(const Float& a, const Float& b)
  60. { return Float(a.num * b.num, a.mes * b.mes); }
  61.  
  62.  
  63. friend int Non_Zero(const Float& f, float i)
  64. { double eps =  i * f.mes * eps0;
  65.  
  66.   if (fabs(f.num) > eps)
  67.      return NON_ZERO;
  68.   else
  69.      if (eps < 1)
  70.         return ZERO;
  71.      else
  72.         return NO_IDEA;
  73.  }
  74.  
  75. friend int Sign(const Float& f)        { return (f.num > 0) ? 1 : -1; }
  76. friend double Lmax(const Float& f)     { return f.mes; }
  77. friend double Eps(const Float& f, int i ) { return i * f.mes * eps0; }
  78.  
  79.  
  80. };
  81.  
  82. inline int Sign(const Int& x)        { return sign(x); }
  83. inline double Lmax(const Int& x)     { return 0; }
  84. inline double Eps(const Int& x, int) { return 0; }
  85. inline bool Non_Zero(const Int& x, float=0) 
  86. { return (Sign(x) != 0) ? NON_ZERO : ZERO; }
  87.  
  88.  
  89. inline int Sign(double x)
  90. { if (x==0) return 0;
  91.   return ( x>0 ? 1 : -1);
  92.  }
  93. inline double Lmax(double)   { return 0; }
  94. inline double Eps(double)    { return 0; }
  95. inline bool Non_Zero(double x, float=0) { return (x != 0) ? NON_ZERO : ZERO; }
  96.